Correlation Effects of Single- Wall Carbon Nanotubes in Weak Coupling 
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Single-wall carbon nanotubes with on-site interaction U <^ t are studied by a controlled renor- 
malization group method. When formulated as a system of Nf flavors of interacting Dirac fermions, 
the effective model at generic filling resembles a set of quasi- ID Hubbard chains. For the undoped 
(A^y, Afj;)-armchair and the {Nx, — A'^i)-zigzag with A^^ — 3n, specific mappings to the 2-chain Hub- 
bard model at half filling with effective interaction u^s = U/Ny and 2U/Nx are found from this more 
general equivalence. The phase diagrams of the (3, 3)-armchair and the (6, -6)-zigzag at generic 
fillings are studied and compared as an example. 



I. INTRODUCTION 

Carbon nanotubes consist of single or mutiple seamless 
shells of graphite sheets with a diameter of about sev- 
eral nanometers and length up to micron scale The 
nanotubes reveal rich varieties of transport properties as 
metals, semi-metals and semiconductors depending on 
their radius and helicity . Their marvelous mechan- 
ical and electrical properties have attracted lots of atten- 
tion in both theory and experiment since their discovery 
in 1991 Recent measurements show that crystalline 
ropes of single-wall nanotubes (SWNT) exhibit metal- 
lic transport properties An individual SWNT has 
been sucessfuUy attached to leads in a two-terminal mea- 
surement and has indicated clear signatures of a ID quan- 
tum wire 1^. The quasi-lD nature of these nanotubes 
necessitates the study of interaction effect, which has 
pronounced influences in low dimensions. A lot of the- 
oretical work focuses on the band structure calculations 
The (A^, Af)-armchair tube with an on-site 
interaction U has been mapped to the 2-chain Hubbard 
model with a reduced interaction strength UcS — U /N 
at or near half filling by L. Balents and M. P. A. Fisher 
1^^ . Similar approaches with more general interactions 
are studied independently by Y. Krotov, D. H. Lee, and 
S. Louie jl^. The correlation effect at generic filling has 
been addressed so far. In this paper, we study a SWNT 
at generic filling with an on-site interaction, which turns 
out similar to the Hubbard chains. 

A SWNT is made by rolling up a 2D graphite sheet 
along the direction of a superlattice vector (TV, M) and 
then identifying all sites related by this vector (for exam- 
ple, see Fig. 1). In this paper, we focus on two specific 
types of SWNT's, the {N, 7V)-armchair and the {N, -N)- 
zigzag. There are equivalent wrappings due to the sym- 
metry of a honeycomb lattice, for example, {N,0) and 
(0,iV) superlattice vectors also produce zigzag tubes. In 
the weak coupling limit U ^ t/N,t±/N, it is natural 
to diagonalize the hoppings and focus on the low-energy 
parts of the whole band structure first. The chemical 
potential cuts through the bands at Nf pairs of Fermi 
points, which enable us to formulate the low-energy the- 
ory as Nf flavors of interacting Dirac fermions. The 
Fermi velocities of these fermions can be computed di- 



rectly from the band structure. Finding all four-fermion 
vertices generated by the on-site interaction proves to be 
trickier. Because of the on-site character, the resulting 
vertices are independent of the underlying lattice struc- 
ture and are shown similar to the Hubbard chains. To 
be precise, interactions of a (iVy, A'j,)-armchair tube at 
the filling, where Nf pairs of Fermi points are present, 
is shown equivalent to the 2iVy-chain Hubbard model 
(with periodic boundary condition) in the Nf partially 
filled band region. Similarly, interactions of a {Nx, —N^)- 
zigzag tube are the same as in the AT^-chain Hubbard 
model except for the absence of the umklapp interactions 
in the ky direction (see below). Knowing the Fermi ve- 
locities and initial couplings of interactions, the phase di- 
agram is determined through the renormalization group 
approach developed previously in Ref . . 

Since the most interesting (feasible) filling is at or near 
half filling, specific mappings for nanotubes to two-chain 
systems based on this more general equivalence are stud- 
ied. We found that the undoped (iVj,, A^j,)-armchair and 
the {Nx, — Nx:)-zigzsLg with A^^, = 3n are equivalent to the 
two-chain Hubbard model at half filling with effective in- 
teraction iteff = U/Ny and 2U /N^ respectively. Unlike 
the previously mentioned 2A^j,- and ATj-chain Hubbard 
model mappings, which only consider interactions, the 
mappings to the two-chain models takes into account in- 
teractions as well as Fermi velocities (band structure). 
In fact, detailed calculations show that couplings of the 
A/-chain Hubbard with an on-site interaction U in the 
Nf = 2 region are identically the same as the 2-chain 
model with reduced coupling strength 2U /N (see below). 
The possibility of finding the mappings to the 2-chain 
model relies crucially on the fact that the armchair and 
the zigzag (A^a, = in) have the same Fermi velocities 
(band structure) as the real 2-chain model at half filling. 

The phase diagrams of the (3,3)-armchair and the (6,- 
6)-zigzag both of which resemble a 6-chain system, are 
analyzed by the RG approach. So far, we haven't ad- 
dressed the point that these tubes in real life are made up 
of "curved" graphite sheets. While this curvature effect 
does not change the low-energy excitations of armchair 
tubes, it creates a minor gap Ac^r ^ t/N"^ in zigzag tubes 
even without interactions (see below) . If the interaction 
U is much smaller than this minor gap, the (6,-6)-zigzag 
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tube is simply an insulator. In the following analysis, we 
focus on the more interesting and non-trivial limit where 
the interaction U is larger than the minor gap caused by 
the curvature effect. Under this assumption, the phase 
diagrams of the armchairs and the zigzags can be com- 
pared at different fillings. At half filling, both of them 
acquire charge and spin gaps caused by interactions and 
thus are expected to be Mott's insulators at zero tem- 
perature instead of metals predicted by simple band cal- 
culation. Upon slight doping into the Nf = 2 region, 
both develop a d-wave quasi-long ranged correlation of 
pair fields. It is a one-dimensional analogy of d-wave su- 
perconductors (SC's). The numerical results of the RG 
analysis indicate that the armchair has a higher critical 
temperature than the zigzag. Larger doping makes the 
armchair (zigzag) enter Nf=6 (4) regime where super- 
conductivity correlations are destroyed. 

The paper is organized as following. In Sec II, nan- 
otubes at or near half filling, i.e. Nf = 2, are studied 
and mapped to equivalent 2-chain models. In Sec III, 
nanotubes at generic fillings are considered. By formu- 
lating into Nf fiavors of interacting Dirac fermions, the 
equivalence of couplings for nanotubes and the Hubbard 
ladders is established. In the last section, implications 
of RG analysis and other theoretical approaches are dis- 
cussed. Long-range interaction effects and open ques- 
tions are also discussed, followed by a brief summary of 
the main conclusions. 



II. NANOTUBES AT OR NEAR HALF FILLING 



Following many authors [^|-|4|, we consider a 2D 
graphite sheet first and then take account of finite ra- 
dius and curvature effects of real tubes. A single sheet 
of graphite is composed of carbon atoms on the sites of 
a honeycomb lattice which is a triangle Bravais lattice 
with two carbon atoms separated by distance d in each 
unit cell as shown in Fig. 1. The sp2 electrons of each 
atom form three bonds with their neighbors and the re- 
maining single Pz electron can tunnel around. Since the 
honeycomb lattice is bipartite, the nearest neighbor hop- 
pings can be viewed as tunneling alternatively between 
two triangle sublattices. A tight-binding model at half 
filling (one electron per site) is introduced to describe 
electron hoppings in a graphite sheet 



r£R,a 



-teL(i")c2a(r + a+ +d) 

-teL(r)c2„(r + a_+d) 
-UcL(r)c2„(r + d) + h.c.|, 



(2.1) 



where R = n^a^ + n_a_ is the lattice vector of a tri- 
angle lattice with bases a± — a(±l/2, •\/3/2). The dis- 
placement vector between these two triangle sublattices 



is d = a(0, — 1/-\/3) and Ci(c|) is the fermion annihila- 
tion (creation) operator on the corresponding sublattices. 
The interaction effect is included by considering on-site 
repulsion between electrons 



Hrj = U 



7i,t(r)n4(r) 



(2.2) 



where riia = cl^c^^ denotes electron density on sublat- 
tice i with spin a. Throughout this paper, only repulsive 
interaction (U > 0) is considered. 




FIG. 1: two dimensional Honeycomb lattice with bases 
a±. Circles and squares represent fermion sites on two 
sublattices separated by a distance d = a/^/3. A (1,1)- 
armchair is obtained by rolling up the sheet from (0,0) 
to (1,1) and end points. Similarly, a (4,-4)-zigzag is ob- 
tained by rolling from (0,0) to (4,-4). 

In the weak coupling regime U <^ t,t±, it is natural to 
diagonalize the hopping Hamiltonian Hq first and study 
the band structure. By going into momentum space, the 
energy spectrum is i?(k) = =F|/i(k)| 

h(k) = 2tcos{k^a/2)e'^y''l^ + t^e-'^y"'^, (2.3) 

where k is the crystal momentum. The ^ signs in the 
energy spectrum denote bonding and antibonding bands. 
Since the undoped graphite sheet is half-filled, the bond- 
ing band is completely filled while the antibonding band 
is empty. One prominent feature of the band structure 
is that the bonding and antibonding bands are separated 
by a finite gap everywhere in the the Brillouin zone (BZ) 
except at two peculiar points K±. They are convention- 
ally called Dirac points because the (gapless) dispersion 
looks relativistic near them, i.e. 



E{\l) =i;|k-K±| 



(2.4) 



For t± = t, the Dirac points are located at K± = 
(±47r/3a) and the Fermi velocity ti = {\/3/2)ta as shown 
in Fig. 2. If t± is slightly different from t, these points get 
shifted a bit along the kx direction. Equivalent points in 
the BZ are related by reciprocal lattice vectors. At half- 
filling (or under slight doping), gapless excitations only 
exist at (or near) the Dirac points and the low-energy 
physics is captured by two flavors of Dirac fermions with 
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the on-site interaction in Eq. 2.2. It is strikingly similar 
to the 2-chain Hubbard model which describes two inter- 
acting Dirac fermions with on-site repulsions. This moti- 
vates us to search for equivalent 2-chain models for nan- 
otubes near half filling. In the following subsections, we 
manipulate the (iVy, A'j,)-armchair and the {N^, —Nj-)- 
zigzag with Nx ~ 3n to show that it is indeed possible 
to map them into equivalent 2-chain models. 



where (±) means, if Ny is an even number, only one sign 
is chosen in order to count the total number of modes 
properly. For notational convenience, we later represent 
the transverse momentum ky by the dimensionless inte- 
ger p. Similarly, since the transverse coordinate y runs 
through Ny values for each fixed x, it is possible to la- 
bel it by a dimensionless integer n. The fermion operator 
Ci{x, n) = Ci{x, y) is then labeled by the integer n defined 
as 



A. Armchair tubes 

The {Ny, A'j^)-armchair has a finite circumference C = 
^/SaNy which leads to the quantization of transverse mo- 
mentum ky around the tube. In addition to the quantized 
momenta, the curvature effect causes the slight differ- 
ence between t, t± which shifts the locations of the Dirac 
points as shown in Fig. 2. Since the quantized momen- 
tum (ky = 0) cuts through the Dirac points, the arm- 
chair without interactions is expected to be a metal. But 
once the interaction effect is taken into consideration, the 
low-energy excitations change dramatically. An armchair 
carbon nanotube with on-site interaction has been shown 
equivalent to the 2-chain Hubbard model at or near half 
filling where only ky = band is gapless ||l2|. Here we 
rederive the same mapping by a more systematic method 
which can also be applied to zigzag tubes later. 



2n 2n 



1 



3 2 10-1-2 



{ 



'(§,0) 



if a; = 



^iy + S^2d) - 2 if x = {m 



(2.6) 



where m,n are integers and n = 1,2, ...,Ny. Since only 
the ky — mode is important for the armchair at or 
near half filling, we want to construct fermion operators 
carrying definite transverse momentum and then discard 
the non-zero momentum components. The required op- 
erator is constructed by partial Fourier transformation in 
the y-direction 



Ci{x,y) — e^^^Ciikxjky) 

^4=5:e'^"^d.(a:,p), (2.7) 



where the normalization constant Ni — N2 — 2NyL/a is 
the total number of fermion sites on each sublattice and 
L is the length of the tube. To obtain the second line in 
Eq. 2.7, we sum over all possible kx and replace y and 



ky by integers n and p respectively. This transformation 
can be written down in a more compact form 



Ci{x,n) = Sil^d^{x,p), 
dz{x,p) = S^yJ^Ci(x,n), 



(2.8) 
(2.9) 



(a) (b) 

FIG. 2: Dirac points in the Brillouin zone. Dark circles 
and squares indicate the locations of the Dirac points for 
t± — t, and grey symbols for t± < t. Dashed lines show 
the allowed quantized momenta in (a) the (3,3)-armchair 
and (b) the (6,-6)-zigzag. The conventional hexagonal 
zone is mapped into rectangle one for (a) the armchair 
and (b) the zigzag. 

To proceed with the mapping, the conventional hexag- 
onal BZ needs to be rearranged into a rectangle. With 
several parts shifted by reciprocal lattice vectors, it is 
done as shown in Fig. 2(a). It is crucial to make such 
a change and the reason will become clear later. With 
the choice of this particular BZ, referred to as BZ(arm) 
later, the allowed transverse momentum takes Ny integer 
values 



_ 2n p 



p = 0,±l,±2,...,(±)[^], (2.5) 



where summation on doubly appeared indices is implied. 
The transformation matrix 5'^*'-* is 



fN, N„ ' 



(2.10) 



This transformation is only possible in the first place 
because the summation over all momenta in the BZ, 
, can be separated into two independent sums, 
X]fe X]fe , which is a direct result of the rectanglar shape. 
Secondly, with the specific choice of the BZ(arm), the mo- 
mentum ky runs through Ny values (given in Eq. 2.5) as 
the transverse coordinate y. As a result, the transforma- 
tion matrix S'^^^ is unitary which preserves the fermionic 
commutation relations between di(x,p) 



{di{x,p),dj{x' ,p')} 



Sij Spp' dxx' 



(2.11) 



Since the Dirac points lie on the ky = Q line, all modes 
with non-zero transverse momentum are gapped and can 
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be integrated out in th e we ak coupling. Keeping the zero 
momentum mode, Eq. 2.8 becomes 



Ci{x,n) 




di{x, 0) = 



1 



= d,;(x), 



(2.12) 



where the dependence of the coordinate y completely 
drops out. Wi thin this approximation, the Hamiltonian 
in Eq. 2.1, 2.1 is mapped to a 2-chain system 



Ho=J2{- td\„ix)d2c.ix - b) 

- td\^{x)d2a{x + b) 

- t±d\^{x)d2a{x) + h.c. L 



(2.13) 
(2.14) 



-^dj(a;) = -0(2;, y), 



(2.16) 



where y = 1, 2 labels different chains. We diagonalize the 
hopping Hamiltonian in the continuous version and ob- 
tain bonding and anti-bonding bands which are denoted 
by ipi{x) with i = ky/n — 0,1. Filling up the bands to the 
chemical potential energy and linearizing the dispersion 
around the Fermi points, the fermion fields in bonding 
and antibonding bands are approximated by two pairs of 
chiral fields as 



-ikp-iX 



(2.17) 



where kpi is the Fermi momentum of each band. The 
effective theory contains two flavors of Dirac fermions 
described by the Hamiltonian 



Wo = ^ Vi 



where b = a/2 is the lattice constant. The hopping term 
may not look familiar at first sight. By depicting hop- 
pings between two sublattices explicitly in Fig. 3(a), it is 
clear that it describes the 2-chain Hubbard model with 
lattice constant b — a/2. It is interesting that the inter- 
action strength reduces to U/Ny. The 1/Ny factor can 
be thought as coming from delocalized electrons around 
the tube since they would then have a probablility 1 /Ny 
to occupy the same site. We emphasize that this effec- 
tive model is only correct when the armchair is at or 
near half filling, such that the other bands with non-zero 
transverse momentum are still gapped. Upon larger dop- 
ing, the interaction effect is more complicated and will be 
studied in the next section. 



t 



t 



(a) (b) 

FIG. 3: Effective 2-chain model for (a) the armchair and 
(b) the zigzag. Squres and circles indicate fermion sites 
on two sublattices. 

The 2-chain Hubbard model has been studied with RG 
techniques by many authors. Here we adapt the approach 
and notations in Ref. ||l^ and briefly review the strategy. 
First of all, the fermions di{x) are relabeled b y cha i n in- 
dices in stead of the sublattice indices in Eqs. p. 13 , ^.14 
Since the coordinate x = mb is an integer multiple of 6, in 
the continuous limit, the sum is replaced by an integral 
and the fermion fields acquire a rescaling factor 



dx 



(2.15) 



(2.18) 

where Vi is the Fermi velocity of each band. There are 
many possible interactions between these Dirac fermions 
but only some of them are relevant to the change of low- 
energy excitations. For example, an interaction with 
right (left) moving fields alone only renormalizes the 
Fermi velocities and can be ignored in the leading or- 
der of U. On the other hand, interactions which couple 
right and left moving fields together can possibly open 
up a gap. To write down these interactions explicitly, 
the SU(2) scalar and vector currents are introduced 



1 i f'^^n'' 



(2.19) 
(2.20) 



where cr is Pauli matrices and e"'^ is the Levi-Civita an- 
tisymmetric tensor with the convention e^^ = — e^^ = 1. 
All possible interactions which are relevant to gap open- 
ings are represented succinctly by these currents 



+ f 7^ 7^ - 



3V 

V V ij 



h.c. 



(2.21) 
(2.22) 



where fij and Cij denote the forward and Cooper scat- 
tering amplitudes respectively, between bands i and j. 
Umklapp interactions vjj fl^ only appear at half filling. 
Summation on i,j is implied. Because of various sym- 
metries, these couplings are symmetrical. Since la = 
due to the anti-symmetric nature of lij, its correspond- 
ing coupling vjf is always set to zero. Expressing the 
the on-site interaction in Eq. 2.14 in terms of t hese Dirac 
fermions, we can compare the results with Eqs. ^.21 , 2.22 . 
The initial couplings generated by the on-site interaction 
are 
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Ub 



Ub 



^^01 - ^^^^ - 4^ 



0, 



(2.23) 
(2.24) 



where i,j = 0, 1 are the band indices. The lattice con- 
stant b appears in the continuous hmit because these cou- 
pUngs are dimensionless. 

With the knowledge of Fermi velocities and initial 
couplings, we employ the RG equations derived in Ref. 

L5| to determine which couplings are (marginally) rel- 
evant and grow into the strong coupling regime when 
short-ranged fluctuations are integrated out. The phase 
diagram is then determined by these relevant couplings 
through the bosonization technique. Numerical results 
show that, at half filling, the relevant couplings gap out 
all spin and charge excitations and lead the system to a 
COSO Mott's insulator, where CnSm means there are n 
gapless charge and m gapless spin modes. The spin and 
charge gaps change the simple band calculation predic- 
tions of the transport properties at low temperature. For 
instance, the resistivity blows up exponentially when the 
temperature is cooled down below the gap scale. Upon 
doping away from half filling, the umklapp interactions 
are no longer allowed. The Cooper scatterings induce a 
net attrative interaction between electron at low-energy 
scale and lead to electron pairings. According to the nu- 
merics, the system goes into the CISO phase with a gap 
function of d-wave symmetry. The spin gap comes from 
the singlet pairing and the only gapless charge mode is 
the usual U(l) phase in a superconductor. 

As one might notice that the results of the RG anal- 
ysis only depends on two factors: the Fermi velocities 
calculated from the band structure and the initial cou- 
plings generated by the on-site interaction. An alterna- 
tive interpretation of the effective model for the armchair 
is available. The Fermi velocities can be calculated di- 
rectly from the band structure of the armchair without 
map ping to any effective models. As to the couplings in 



Eqs. 2.23, 2.24 , they are similar to the iV-chain Hub- 



bard model with lattice constant b, (e.g. cfj — 2Ub/N). 
By comparison, the initial couplings for the (Ny,Ny)- 
armchair in the Nf — 2 region are identical to those of 
the 2A'j,-chain Hubbard model with two partially filled 
bands, where Nf is the number of Dirac fermions in the 
nanotubes. It just happens that the band structure for 
fcy = in the armchair has the same dispersion as in a 
real 2-chain system, which enables us to find the map- 
ping described before. This new interpretation has the 
advantage that the equivalence of the couplings holds at 
generic fillings, where an effective 2-chain model is im- 
possible. The connection between the {Ny, 7Vj,)-armchair 
and the 2iVy-chain Hubbard model is elucidated in detail 
in the next section. 



B. Zigzag tubes 

For an {N^, —Nx)-zigza.g tube, the circumference is 
C = Nxtt and the momentum kx is quantized. The 
band structure shows interesting results even when in- 
teractions are negligible. For now, we ignore the cur- 
vature effect and set t± = t. For = 3n, the Dirac 
points lie on the allowed momenta (see Fig. 2(b)). The 
system has gapless excitations and thus is a metal. For 
Nx ^ 3n, the allowed momenta do not cut through the 
Dirac points and the (filled) bonding and (empty) anti- 
bonding bands are separated by a finite gap. Thus, the 
system is a semi-conductor. Due to curvature effects, 
the difference between transverse and longitudinal hop- 
pings is of order I /N"^ 7 and it causes a gap A also 
of order l/iVj even when A^^, = 3n. However, we find 
later that the effective interaction is of order l/N^. It is 
then reasonable to ignore the minor gap caused by curva- 
ture and treat (3n, — 3n)-zigzag tubes as metallic. Since 
Nx 3n tubes are already insulating without interac- 
tions, the inclusion of interaction effect does not change 
the low-energy physics dramatically. We focus on the 
more interesting case of A^^; — 3n and explore whether 
the interactions can indeed cause gaps. 

For the (3n, — 3n)-zigzag, a mapping to a 2-chain 
model, similar to the armchair, is found by partial Fourier 
transform in x. To perform the transformation, as for 
the armchair, a rectanglar BZ is again necessary. How- 
ever, the BZ(arm) is not suitable because it contains 2Nx 
quantized kx which is twice the number of modes neces- 
sary to make the transformation unitary. Another BZ 
is chosen as shown in Fig. 2(b), which is referred to as 
BZ(zz) later. With this choice, the momentum kx has 
Nx allowed values: 



kx — 



27r p 



N 

p = 0,±l,±2,...,(±)[-^]. (2.25) 



For each fixed coordinate y, the transverse coordinate 
X runs through Nx possible values. The fermion field, 
c(n, y) = c{x, y), is labeled by the integer n defined as 



x/a 
x/a - 



if y 
if y 



m\/3a, 

(m + |)\/3a. 



(2.26) 



where m,n are integers and n — l,2,...,Nx. The 
fermion field di (p, y) which carries definite momentum 
kx is constructed by partial Fourier transformation in x 



Ci{n,y) = Sl^jld^{p, y), 
Mp^v) = Sln^c^(n,y). 
The transformation matrix S^^^ is 



-;= explj— — ao). 

fWx ^ ^ Nx ' 



(2.27) 
(2.28) 

(2.29) 



Since S^^^ is unitary, the fermion field di(p,y) obeys the 
canonical anti-commutators. For tj_ = i, the Dirac points 
lie on two specific momenta, kx — ±27r/3a, and other 
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bands acquire major gaps. Even with the inclusion of 
curvature effects, the gaps of the = ±27r/3a bands is 
only of order which is smaller than those of the 

other bands. In the weak coupling, it is legitimate to 
integrate out all other gapped modes and Ci{x,y) only 
contains components near the Dirac points 



1 



q=± 



(y)e 



(2.30) 



where d±i{y) = di{±Nx/iTy). With in t his approxima- 
tion, the hopping Hamiltonian in Eq. 2.1 is simplified to 



two chains without interchain hoppings 

i/o= E {-tdlAy)d,2{y + b' + 5) 

y,g=± ^ 

-t±dl^{y)dg2{y - b' + S) +h.c. 



(2.31) 



where q — ± is the chain indices. The lattice constant is 
b' = \/3a/4 and the distortion distance is (5 = ^/3a/12. 
The distortion S can be gauged away by shifting one of 
the sublattices. Performing a gauge transformation to 
shift the second sublattice by 6, we get 



DU,iiy)D^d,,iy), 
DU,2iy + 6)D^dM. 



(2.32) 
(2.33) 



In momentum space, this gauge transformation only adds 
a phase to dq2{k) dq2{k)e^''^ and does not change the 
band structure at all. As to the interaction, since it is on- 
site, it is also invariant under this gauge transformation. 
The hopping Hamiltonian becomes 

Ho = Y.{-tiy)dliy)d^iy + h') + h.c.}, (2.34) 



with the hopping t{y) alternating along the y-direction 
as 

t{y) = \{{t + t^) + {-i)y{t^t^)]. (2.35) 

For t±_ 7^ t, the hopping Hamiltonian represents a dimer- 
ized chain as shown in Fig. 3(a) and has a gap A ~ 
(t^ — t) while, for t±_ = t, the gap vanishes. This is in 
agreement with our previous conclusion for the zigzag. 
As argued before, the minor difference of order be- 
tween and t± caused by curvature effects is negligible 
because the interactions are order of 1/Nx- Thus we set 
i = in the following. It takes a little algebra to write 
down the interactions between these two chains and ob- 
tain 



"9T (y)"9i (y) + '^-jt iy)nqi {y) 



dl^'iy)dq^iy)dli^y)dqi{y) 



(2.36) 



where the notation a = —a is adapted through out the 
paper. It is interesting to compare the effective 2-chain 
models for the armchair and the zigzag. For the arm- 
chair, the on-site interaction remains simple but there 
exist hoppings between two chains, while, for the zigzag, 
there is no hoppings between chains but the interactions 
are complicated. The effective model obtained here does 
not look like the 2-chain Hubbard model yet. Follow- 
ing the method for the armchair, we take the continuous 
limit and linearize the spectrum near the Fermi points. 
The fermion fields in these two chains are approximated 
by two pairs of chiral fermions as 



My)-^Rqiyy'''''' + ^LMe 



KFqV 



(2.37) 



Once more we have an interacting theory with two fla- 
vors of Dirac fermions. Because of the absence of in- 
terchain hoppings, these two bands are degenerate and, 
thus, always have the same Fermi velocities. The al- 
lowed interactions are the same as those of the armchairs 



in Eqs. 2.21, 2.22. The initial couplings determined from 
the interactions in Eq. 2.36 are 







4rf. = 


Ip 




Ub' 
2Nx 



fij 



2Ub' 

Nx ' 



0. 



(2.38) 
(2.39) 



The couplings are exactly the same as for the 2-chain 
Hubbard model with reduced interaction u^g = U/N^ 
and lattice constant b' = \/3a/4. At half filhng, the 
Fermi velocities of Dirac fermions for the armchair and 
the (3n, , — 3n)-zigzag are identically the same. As a re- 
sult, the effective theory of the zigzag is again described 
by the 2-chain Hubbard model with reduced interaction 
U/Nx and lattice constant h' — ■\/3a/4. Upon doping, 
the Fermi velocities of the zigzag remain degenerate and 
start to deviate from the Hubbard model. However, nu- 
merical results show that the phase is still CISO as in 
a doped armchair tube. It is interesting to compare a 
(3iV, 3iV)-armchair and (6A^, — 6A^)-zigzag because the ef- 
fective interaction strength Ucs ~ U/3N is the same for 
both. However, since the lattice constant of the armchair 
b is slightly larger, the initial couplings of the armchair is 
a bit larger than those of the zigzag. By integrating the 
RG flow numerically, the results imply that the critical 
temperature of the armchair is higher than the zigzag due 
to the difference in the initial couplings from difference 
lattice constants. 



III. NANOTUBES AT GENERIC FILLINGS 

So far, we studied the carbon nanotube at fillings de- 
scribed by two flavors of interacting Dirac fermions and 
found it is closely related to the 2-chain Hubbard model. 
Upon larger doping, the chemical potential cuts through 
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more pairs of Fermi points and, in general, there are Nf 
flavors (more than two) of Dirac fermions participating in 
the interactions. For example, the (3, 3)-armchair enters 
Nf = 6 region upon larger doping as indicated in Fig. 4, 
while the (6, — 6)-zigzag enters Nf = 4 region as shown in 
Fig. 5. The strategy used to analyze them with RG tech- 
niques involves two steps: calculate the Fermi velocities 
from the band structure, and figure out all possible inter- 
actions with the initial couplings generated by the on-site 
interaction. The first step is easy because the band struc- 
ture of nanotubes is already obtained in Eq. 2.3. The 



second step is more challenging. We are going to show 
that the interactions between these Nf flavors of Dirac 
fermions are almost the same as in the A'^-chain Hubbard 
model in the Nf partially filled band regime. In the fol- 
lowing, we focus on the (3, 3)-armchair in the Nf = 6 
region and the (6, — 6)-zigzag in the Nf — 4 region as ex- 
amples. The generality of this equivalence is transparent 
from these two typical illustrations. 





(a) (b) 

FIG. 4: The (3,3)-armchair in iV/ = 6 partially filled 
band region. The band structure near E' = is shown 
in (a) and the corresponding Fermi points are also plot- 
ted in the BZ in (b). The thick line represents the two 
degenerate bands with transverse momentum p = ±1. 
For each quantized momentum p, there are two flavors 
of Dirac fermions participating in the interactions. The 
solid (empty) arrow lines indicate the umklapp coupling 
^12 (""12) addition to the Cooper and forward cou- 
plings. 

We study the armchair flrst. For convenience, the 
sublattice index i — 1,2 in Ci{x,y) is viewed as the z- 
coordinate. The corresponding momentum p' = k^/n = 
0, 1 labels the bonding and anti-bonding bands respec- 
tively. Perform a partial Fourier transformation in 
BZ(arm) for both y, z to diagonalize the hoppings around 
the tube, we get 



Ci{x,n) — S^JS^}cp{x 



(3.1) 



where p = {p,p') is the momentum in the y, z directions. 
The transfer mat rix S'^^^ is defined analogous to S*^^^ de- 
fined in Eq. |2T^ 



S. 



(^) 



:exp(i— a6), 



(3.2) 



where the normalization constant Nz = 2 accounts for 
two possible values of z{kz). In the continuous limit, the 
fermion field is rescaled as 



Cp(a;) = Vbij)p{x). 



(3.3) 



We linearize the dispersions near each Fermi point and 
describe the system by Nf = 6 fiavors of Dirac fermions. 
Since the system is doped, only the bonding band {p' — 
0) participates in the interactions. Since the anti-bonding 
band (pi' = 1) never appear in the following analysis, the 
p' index is dropped for notational simplicity. To make 
the equivalence to the Hubbard model transparent, the 
fermions are labeled in a particular way: 



+ Vflp. (x)e*'=^-*" + VL„-„.(a^)e 



-ikj; 



(3.4) 



where p* = p — sgn{p)Ny. The labeling of these fermions 
is given explicitly in the Fig. 4. It is straightforward 
to check that the Cooper and forward scatterings are al- 
lowed between any two of the fermions. However, since 
there are two flavors of Dirac fermions for each ky^ ad- 
ditional interactions are also allowed. Two examples are 
shown in Fig 4(b). These are the umklapp interactions 
in the ky direction 



-yiM _ Ip jR jL _ la jR . jL 



2p jR jL _ -la jH . jL 



,2cr jR 



(3.5) 



The name comes from the fact that if the BZ(arm) is 
shifted back to the the conventional one, the transfer 
transverse momentum associated with each vertex equals 
the reciprocal lattice vector An/^/Sa. T he al l owe d inter- 
actions are therefore described by Eqs. 2.21 , 3.5 as in a 



doped (even-chain) Hubbard model with periodic bound- 
ary conditions. We need to compute the initial cou- 
plings generated by the on-site interaction to complete 
the equivalence. The on-site interaction is expressed in 
terms of the Dirac feermions as 



Hu = U^ n^{x,y,z)ni{x,y,z) 

x,y,z 



323334, 



(3.6) 



where Iq is the momentum index of ky and run through 
Nf — 6 values. The ja labels the momentum and is 
always zero because only the boning band is interactiing. 
The chirality of the fermions are denoted by Pi^ = R,L. 
The phase associated with each vertex is the transfer mo- 
mentum in x-direction given by. 



-Piikpii 



Pi^kpis + Pi^kpu- (3.7) 



Conservation of crystal momentum requires Qx equal to 
one of the reciprocal lattice vectors. If the momentum is 
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not conserved, the phase associated with the vertex os- 
cihates and the vertex is irrelevant in the leading order 
analysis. The summations of the transformation matrices 
in the y, z directions are 



id) ^ 



id)* aid) aid)* aid) _ 



QdiGd 



(3.8) 



where d = y,z and Qd — (— ii + ^2 — ^3 + ii) is the trans- 
fer momentum and Gd = '2nNd is the reciprocal lattice 
vector. The J-function enforces the crystal momentum 
conservation requirement. Putting everything together, 
the interacting Hamiltonian, in the continuous limit, is 

Hu = (^) j rfa;V'l..^,,tV'P,,»,T^k»3i^PM'4i' (3-9) 

where summation over all possible vertices Pi^,ia, con- 
strainted by momentum conservation, is implied. These 
allowed vertices which can be cast into current-current 

have been discussed 



interactions in Eqs. p. 21 
fore 



3.5 



be- 



The key point here is that the prefactor Ub/2Ny 
obtained above is indeed exactly the same as that of the 
2Ai"y-chain Hubbard model with lattice constant b — a/2 
in the Nf filled band region. Therefore the equivalence 
of interactions between the (iVy, A^y)-armchair and the 
2A^j,-chain Hubbard model becomes clear. 

Knowing the above equivalence, the (3,3)-armchair in 
the Nf = 6 region can be studied by integrating the 
RG equations numerically. Two of the Dirac fermions 
at fcy = become completely decoupled and contribute 
two gapless charge and two gapless spin modes (C2S2). 
One recalls that these two Dirac fermions form d-wave 
pairs and result in a CI SO superconducting phase in the 
Nf = 2 region. Relevant interactions between the re- 
maining four flavors of fermions are too complicated to 
be bosonized. Therefore, the phase of the armchair in 
the Nf = 6 region remains an open question for future 
study. 



Now we turn to the zigzag. Perform a partial Fourier 
transformation in both a;, z to diagonalize the hoppings 
around the tube: 



c{y,m',n') = s';^}^si^,lcp{y) , 



(3.10) 



Since the coordinate y — 2nb' is an integer multiple of 
26', in the continuous limit, the fermion field is rescaled 
as 



Cp(y) = V2b'^Pp{y). 



(3.11) 



Linearizing the spectrum near the Fermi points results in 
following four flavors of Dirac fermions 



— ikFpX 



(3.12) 



The Cooper and forward scatterings are again allowed 
because the transfer momenta of these vertices are con- 
served. The ky umklapp interactions are not allowed be- 
cause the Fermi momenta do not sum up to the reciprocal 
lattice vector as shown in Fig. 5. The initial couplings 
of these interactions are 



H, 



,Ub\ 



This interacting Hamiltonian is the same as the Nx-c\i&in 
Hubbard model with lattice constant b' in the Nf par- 
tially filled band region, except that the ky umklapp in- 
teractions are missing. 

With this equivalence, the (6,-6)-zigzag may be stud- 
ied as follows. Upon doping into the Nf = 4 region, the 
superconductivity vanishes. The bands near the Dirac 
points are decoupled and contribute 2 gapless charge and 
2 spin excitations (C2S2). The remaining two bands form 
the usual "d-wave" pairing and contribute CISO. The fi- 
nal phase then is C3S2 and the system no longer has a 
spin gap. This result shows that the superconductivity 
only exists inside a narrow window of doping for these 
nanotubes. 




(a) (b) 

FIG. 5: The (6,-6)-zigzag in TV/ = 4 partially filled band 
region. The band structure near = is shown in (a) 
and the corresponding Fermi points are also plotted in 
the BZ in (b). The interaction between Band ±1,±2, 
as indicated by solid lines, is not allowed because the 
momentum transfer is not a reciprocal lattice vector. 



IV. DISCUSSIONS AND CONCLUSIONS 

Experimental studies of multi-wall nanotubes show 
that the electrical properties vary from tube to tube, and 
are dominated by disorder and weak localization in stead 
of the tubule structure 0,11^ . The single- wall nanotubes 
(SWNT's) on the other hand have simpler structures and 
are expected to be more helpful to clear up various inter- 
action effects. Recently, they have been synthesized with 
high yields and structure uniformity 0. The (10,10)- 
armchairs dominate these bulk samples |[19|] . Tubes dur- 
ing the production process form a 2D triangular lattice 
(crystalline ropes) and eventually become mats (three 
dimensional sample with entangled ropes). The (10,10)- 
armchair is predicted to be metallic (within band calcu- 
lation) in agreement with some of the experimental ev- 
idence. Some of the transport measurements show that 
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the resistivity of the ropes and mats increases hnearly 
with temperature in the high temperature regime, while 
at low temperature, one observes a crossover to negative 
dp/dT 0. However, the experimental results are not yet 
conclusive. From purely electronic (Coulomb) interac- 
tions, it was shown that the umklapp interactions at half 
filling cause a linear increase of resistivity with temper- 
ature in the region above the Mott insulating gap, and 
an exponential increase upon cooling below the gap . 
Another approach explains the linear resistivity at high 
temperature by the study of thermal shape fluctuations 
of these tubes [0|o|. These two approaches produce 
the same predictions at high temperature. However, at 
low temperature, the resistivity caused by phonons still 
decreases linearly with temperature in contrast to the ex- 
ponential increase due to umklapp interactions. A pre- 
cise resistivity measurement at low temperature is nec- 
essary to distinguish these two contributions. However, 
the observed dependence of this crossover on the samples' 
morphology and quality suggests that disorder or (and) 
couplings between single tubes may dominate the low 
temperature behavior. Because the present experimental 
results are not yet clear in this regime, it will be impor- 
tant to carry out more experiments to understand the 
origin of this turnover in the resistivity measurements. 



Upon doping, both of the {N, A^)-armchair and 
{3N, — 3A^)-zigzag with on-site interactions become su- 
perconductors (SC's) as discussed in the text. Recent 
studies show that this SC phase is stable with inclusion 
of moderate nearest neighbor interactions pl| ]. It im- 
plies that the SC phase probably always exists for generic 
short-range interactions. So far, there is no theoretical 
work which includes the long-range interaction effect. It 
would be interesting to see how the unscreened Coulomb 
repulsion affects this superconducting behavior in the fu- 
ture study. 



In conclusion, we represent nanotubes at generic fill- 
ings as Nf flavors of interacting Dirac fermions whose ve- 
locities are determined by the band structures. With spe- 
cific on-site type interactions for the (A'j,, iVj,)-armchair 
and the {N^, — iVa;)-zigzag, we found the vertices equiva- 
lent to those of the 2Ny- and A'a;-chain Hubbard model 
respectively, which can be analyzed by the RG method in 
the weak coupling. With this equivalence, mappings from 
nanotubes to the effective 2-chain model with reduced 
couplings are achieved at or near half filling. The (3,3)- 
armchair and the (6,-6)-zigzag are studied and compared 
as an example because both of them resemble 6-chain sys- 
tems. At half filling, they are both Mott's insulators and 
become superconductors upon slight doping. A larger 
doping to Nf > 2 regime destroys the superconductivity. 
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